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We develop a theory of electron transport through quantum dots that are weakly coupled to 
ferromagnetic leads. The theory covers both the linear and nonlinear transport regime, takes non- 
coUinear magnetization of the leads into account, and allows for an externally applied magnetic field. 
We derive generalized rate equations for the dot's occupation and accumulated spin and discuss the 
infiuence of the dot's spin on the transmission. A negative differential conductance and a nontrivial 
dependence of the conductance on the angle between the lead magnetizations are predicted. 
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I. INTRODUCTION 

The study of spin-dependent transport phenomena has 
recently attracted much interest and led to new device 
functionalities based on the manipulation of the spin de- 
gree of freedomi^ A prominent example, which has al- 
ready proven technological relevance, is the spin valve 
based on either the giant magnetoresistance effecls^i in 
magnetic multilayers or the tunnel magnetoresistance^ in 
magnetic tunnel junctions. In the latter case, the tunnel- 
ing current between two ferromagnets depends on the rel- 
ative orientation of the leads' magnetizations. The maxi- 
mal and minimal transmission is achieved for the parallel 
and antiparallel configurations, respectively. At interme- 
diate angles </> of the relative magnetization direction, 
the transmission is interpolated by a cosine law. This 
angular dependence, which has been beautifully demon- 
strated experimentally,-^ can be easily understood within 
a noninteracting-electron pictured It simply reflects the 
(/)-dependent overlap of the spinor part of the majority- 
spin-electron wave functions in the source and drain elec- 
trodes. 

On the other hand, tunneling transport through nanos- 
tructures, such as semiconductor quantum dots or metal- 
lic single-electron transistors, is strongly affected by 
Coulomb interaction, and a noninteracting-electron pic- 
ture is no longer applicable. Recent works on spin- 
dependent tunneling through these devices include stud- 
ies on metallic islands,^"- granular systems)^ and car- 
bon nanotubes^'' coupled to ferromagnetic leads as well 
as spin-polarized transport from ferromagnets through 
quantum dotsi"i^^i"i^^i^^i^^i^^i^^i^Q The focus of this pa- 
per is on quantum-dot spin valves, i.e., quantum dots 
attached to ferromagnetic leads, with non-collinear mag- 
netizations, see Fig. Q] 

The much richer transport behavior of quantum-dot 
spin valves, as compared to single magnetic tunnel 
junctions, relies on the possibility to generate a non- 
equilibrium spin accumulation on the quantum dot, de- 
pending on system parameters such as gate and bias volt- 
age, charging energy, asymmetry of the tunnel couplings, 



and external magnetic field. This opens the potential 
of a controlled manipulation of the quantum-dot spin, 
detectable in transport. To some degree, there is a rela- 
tion between the system under consideration and a single 
magnetic-atom spin on a scanning tunneling microscope 
tip. For the latter, precession of the single spin in an 
external magnetic field has been detected in the power 
spectrum of the tunneling current^^*^ 




FIG. 1: Quantum-dot spin valve: a quantum dot is weakly 
coupled to ferromagnetic leads. The magnetization directions 
of the leads enclose an angle <j>. 

In a recent Letter by two of us,— we predicted that 
the interplay of spin-dependent tunneling and Coulomb 
interaction in quantum-dot spin valves gives rise to an 
interaction-driven spin precession, describable in terms 
of an exchange field and detectable in a nontrivial (f> 
dependence of the linear conductance. The aim of the 
present paper is to extend this discussion to a complete 
theory of transport through quantum-dot spin valves in 
the limit of weak dot-lead coupling;^: The theory covers 
both the regimes of linear and nonlinear response, and 
is easily generalized to include an external magnetic field 
and/or spin relaxation processes within the quantum dot. 
The quantities of central interest are the magnitude and 
direction of the spin accumulated on the quantum dot. 
They can be determined by generalized kinetic equations 
within a diagrammatic real-time transport theory for the 
dot's density matrix. To obtain a more intuitive and eas- 
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ily accessible framework for the dynamics of the dot's 
occupation and spin, however, we transform the kinetic 
equations into rate equations that exhibit similarities to 
Bloch equations of a spin in a magnetic field. The deriva- 
tion of these Bloch-like rate equations is the first central 
result of this paper, presented in Sec. IIIII after defining 
the model in Sec. ^1 They provide a convenient starting 
point for the discussion of how the dot state depends on 
the system parameters. The second major achievement 
of this paper is to relate the current through the device 
to the dot's occupation and average spin, which is done 
in Sec. 

An analysis of the obtained results is presented in 
Sec. 13 We find that the quantum-dot spin, and thus 
the current, behaves quite differently in the linear- and 
nonlinear-response regime. In linear response, the above- 
mentioned effective exchange field rotates the quantum- 
dot spin such that the spin-valve effect is reduced. This 
contrasts with the nonlinear-response regime, in which 
the dot spin tends to align antiparallel to the drain elec- 
trode, leading to a spin blockade. Together with the ef- 
fective exchange field this leads to a very pronounced 
negative differential conductance. 

We, then, conclude by generalizing our theory to in- 
clude an external magnetic field fSec lVIjl or internal spin 
relaxation in the dot fSec. IVIlll . and by interpreting the 
effective exchange field within the recently introduced^^ 
language of the spin- mixing conductance fSec. IVIlH) . 

II. MODEL 

We consider a quantum dot with a level spacing exceed- 
ing thermal broadening, intrinsic linewidth, and charging 
energy, i.e., only a single level with energy e, measured 
relative to the Fermi energy of the leads, contributes to 
transport. The quantum dot is tunnel- coupled to ferro- 
magnetic leads. We model this system with the Hamil- 
tonian 

H = H,i,t + Hi^ + Hn + HT,L + HT^n. (2.1) 

The first part i?dot — X^o- £cJ.Co. + U n^nj^ describes an 
atomiclike spin-degenerate energy level on the dot plus 
the charging energy U for double occupancy. The left and 
right ferromagnetic leads, r = L, R, are treated as a reser- 
voir of itinerant electrons, Hr = J2ka ^ka^lka'^rka' where 
a — + labels majority and a = — minority-spin electrons. 
Here, the quantization axis for the electron spins in reser- 
voir r is chosen along its magnetization direction n^. In 
the spirit of a Stoner model of ferromagnetism, we assume 
a strong spin asymmetry in the density of states Pr,±(^) 
for majority (-I-) and minority (— ) spins. (To be more 
precise, without loss of generality we assume the direction 
of magnetization to be parallel to the direction of major- 
ity spins.) In the following, the densities of states are 
approximated to be energy independent Pr,±(^) = Pr,±- 
(Real ferromagnets have a structured density of states. 
This will modify details of our results but not change the 



general physical picture.) The asymmetry in the density 
of states is characterized by the degree of spin polariza- 
tion Pr = {pr+ — Pr-)/{Pr+ + Pr~) with < Pr < 1, 

where Pr — Q corresponds to a nonmagnetic lead, and 
Pr = 1 describes a half-metallic lead, carrying majority 
spins only. The magnetization directions of the leads can 
differ from each other, enclosing an angle 4> = <(nL, Ar). 
An applied bias voltage is accounted for in different elec- 
trochemical potentials for the left and right leads. 

Tunneling between leads and dot is described by tun- 
nel Hamiltonians -ffx.L and H^.r- Their explicit form 
depends on the choice of the spin quantization axis for 
the dot states. We find it convenient to choose neither hl 
nor but instead to quantize the dot spin cr J, along 
the z-direction of the coordinate system in which the ba- 
sis vectors ex, and are along hl -I- hr, hl — nR, 
and fiR X fiL, respectively (see also Fig. |21l. The tunnel 




FIG. 2: (a) Scheme of the quantum-dot energy scales. The 
difference of the electrochemical potentials of the left and 
right leads describes a symmetrically applied bias voltage, 
(b) The coordinate system we choose. The magnetization 
directions hl and hr enclose an angle 0. 

Hamiltonian for the left tunneling barrier then reads 

+ aL_(-e+^^/%+e-^/%)+H.c. (2.2) 

and for the right barrier, i?T,R, the same with the re- 
placements L — > R and (j) ~^ ~4'- Here, ti, and iR 
are the tunnel matrix elements. The tunneling rate for 
electrons from lead r with spin (±) is then Tr^±/h — 
2TT\tr\'^ Pr^±/h. Therefore the electronic states in the 
quantum dot with spin (||) acquire a finite line width 

rr = E.=± 2^|ijV2|2p,, , = Y.a=± rr,./2. 

We note that while individual parts of the tunnel 
Hamiltonians do not conserve spin separately, the sum 
of all parts strictly does. Due to the special choice of the 
quantization axis, the lead electrons are coupled equally 
strong to up- and down-spin states in the quantum dot. 
There are, however, phase factors involved, similar to 
Aharonov-Bohm phases in multiply connected geome- 
tries. The formal similarity of the quantum-dot spin 
valve to a two-dot Aharonov-Bohm interferometer^^ is 
visualized in Fig. |21 
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FIG. 3: Graphical representation of the phase factors in the 
tunnel Hamiltonian, Eq. 12.211 . In this representation, the 
quantum-dot spin valve bears similarities to an Aharonov- 
Bohm setup, with the different arms modeling the two dif- 
ferent spin channels t and j, and the Aharonov-Bohm flux 
corresponding to the relative angle of the leads' magnetiza- 
tions. 



III. OCCUPATION AND SPIN OF THE DOT 

In this section we derive rate equations governing the 
state of the dot, i.e., its occupation as well as the mag- 
nitude and direction of an accumulated spin. For this, 
we make use of a diagrammatic real-time approach^'^ to 
describe the time evolution of the reduced density matrix 
for the dot degree of freedoms, which is obtained after in- 
tegrating out the lead electrons. The resulting equations 
for the density matrix elements can be reformulated into 
a more intuitive language of rate equations for the dot's 
occupation and spin. 



A. Reduced density matrix for the quantum dot 

The dynamics of the system can be described by the 
time evolution of the total density matrix, which accounts 
for all degrees of freedom in the dot and the leads. While 
the leads are treated as reservoirs of noninteracting elec- 
trons, it is the dynamics of the dot degrees of freedom 
that determine the transport behavior. This motivates 
the idea to integrate out the leads and arrive at a reduced 
density matrix for the dot degrees of freedom only. De- 
tails about the procedure of integrating out the leads are 
given in Ref . 

The basis of the Hilbert space is given by the states 
X = (empty dot), t (dot occupied with one spin-up 
electron) , | (the same with a spin-down electron) , and d 
(double occupancy of the dot). Therefore, the reduced 
density matrix p^ot is a 4 x 4 matrix with matrix elements 
We get 



-Px2 ~ (^dot)xiX: 



^0° 



Pdot — 





Ft 

pi 







Pi 











p<lj 



(3.1) 



The diagonal, real entries P^ = P^ are nothing but the 
probabilities to find the dot in the corresponding state 
X = 0, t: ii d. The zeros in Eq. (|3.1|) in the off diagonals 
are a consequence of the total-particlc-number conserva- 
tion. The nonvanishing complex off-diagonal elements 
Pj = (Pj)* reflect the fact that an accumulated dot spin 
is not restricted to the z direction but can also have a 
finite X and y component. The density matrix has, thus, 
six real parameters (four diagonal elements and the real 
and imaginary parts of Pj). Due to the normalization 
condition Tr pdot = 1, only five of them are independent. 

The time evolution of the reduced density matrix 
is given by a generalized master equation in Liouville 
spaced, 



--{E^,^E^,)P^l{t) 



dt' 



Xi X2 



P'^Ht') 



(3.2) 



Here, E 
case Eq = 0, E 



. are the energies of the decoupled system (in our 



nels Yf'] it' , t) act as generalized transition rates in Li- 

Xi X2 ^ ' ' *= 

ouville space. They are defined as irreducible self-energy 
parts of the dot propagator on a Keldysh contour (see 
appendix lA 1(1 , and can be expanded in powers of the 
dot-lead coupling strength F. In the limit of interest for 
this paper, i.e., weak dot-lead couphng, only the terms 
linear in F are retained. 

The Markovian approximation is established by re- 
placing P^^^it') in Eq. (E3 with P^/(t). In this case, 
the time integral can be comprised in the definition 
j:^'} XI ^ -f^ fO ^^,j.xi XI QY ^ ^^Yl only discuss re- 

X2X2 J-00 X2X2^ ' ' 

suits for the steady-state limit, i.e., all matrix elements 
are time independent and the above-mentioned replace- 
ment is not an approximation but valid in general. In 
order to emphasize the physical origin of the rate equa- 
tions, however, we will keep all time derivatives although 
their numerical value is zero for our discussion. 



£, and Ed^2e + U). The ker- 



B. Master and Bloch-like equations 

For a more intuitive access to the average dot state, 
we connect the matrix elements of the reduced density 
matrix to the vector of the average spin hS with 5' — 
(5x, Sy, S^) by 



P + Pi 



Sy = t- 



P, 



pi 



(3.3) 



We emphasize that the S describes a quantum-statistical 
average rather than the coherent state of a single spin. 
As a consequence, the magnitude \S\ can adopt any num- 
ber between and 1/2. The average charge occupation 
of the dot is determined by the probabilities Pq, Pi, and 
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Pd, where Pi = P| + P|. The dot state is, thus, char- 
acterized by the set of parameters (Pq, Pi,Pd, S'x, Sy, 6*2), 
in which only five of them are independent, since the nor- 
mahzation of the probabihties requires Pq + Pi + Pd = 1 • 
We evaluate the irreducible self-energies appearing in 



Eq. H3.2(l up to first order in F (see Appendix to get 
a set of six master equations, three for the occupation 
probabilities and three for the average spin. The first 
three are given by 





2prTr 



-/-(£)-/+(£ + [/) 2f-{e + U) 
f+{e + U) -2f-{e + U) 

-f-{e) + f+{e + U) 




(3.4) 



Here, we have used the following definitions. The prob- 
ability to find an electronic state at energy uj in the 
left (right) lead occupied is given by the Fermi function 
/l(r) i^) by (w) = 1 - (w) to find it empty. 

In the case of nonmagnetic leads, Pr = 0, only the 
first line in Eq. H3.4|l survives, and we recover the usual 
master equations for sequential tunneling through a sin- 
gle leveli^ For spin-polarized leads, however, tunneling 
into or out of the dot becomes spin dependent. As a 
consequence, the charge occupation, in general, depends 



on the spin S accumulated on the dot, as indicated by 
the second line of Eq. H3.4|l . The three remaining master 
equations describe the time evolution of the average spin, 

dS fdS\ fdS\ fdS\ 

with 



dS_ 
'dt 

dS_ 
It 

dS_ 

'dt 



rcl 



E 



Pr 



= ^ Br, 

r=L,R 



(3.6) 
(3.7) 
(3.8) 



where we used the definition 

■Kfi J \uj — e — U UJ — e J 

for r = L, R, and the prime at the integral symbolizes 
Cauchy's principal value. We see that three different con- 
tributions lead to a change of the average dot spin. The 
first one, Eq. H3.6|l . depends on the occupation probabili- 
ties only and describes noncquilibrium spin accumulation 
via tunneling to and from spin-polarized leads. This is 
the source term, which is responsible for building up a 
spin polarization in the quantum dot. The second term 
(proportional to S) has the opposite effect. It accounts 
for the decay of the dot spin by tunneling out of the elec- 
tron with given spin or by tunneling in of a second elec- 



tron with opposite spin forming a spin singlet on the dot. 
These tunneling mechanisms yield an intrinsic spin relax- 
ation rate = Er^^^i'^r/h) [f- {e) + f+{e + U)]. 

From Eq. H3.5|l we can derive the magnitude of the spin 
on the dot via 

dJS\_S_ fdS\ _\s\ 

dt -\S\ UJacc ■ ^'-''^ 

The two contributions, Eqs. (|3.6|) and H3.7|) . capture 
the spin dynamics due to spin-dependent tunneling pro- 
cesses transferring electrons into and out of the quantum 
dot. But above equation is also sensitive to the direction 
of S, which depends directly on the third contribution, 
Eq. H3.8|) . of the Bloch-like equation. This term leads 
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to a precession of the dot spin about Bl + -Br, where 
(Bl + -Br) /7 can be viewed as a fictitious magnetic field, 
with 7 — —guB being the gyromagnetic ratio. This ficti- 
tious field vanishes in the absence of Coulomb interaction, 
U = 0, i.e., its origin is a many-body interaction effect. 
Its nature can be conceived as an exchange field due to 
virtual particle exchange with the spin-polarized leads. 
Note that these virtual exchange processes do not change 
the charge of the dot, in contrast to the spin-dependent 
tunneling events responsible for the first two contribu- 
tions to Eq. 1)3. 5|l . Nevertheless, this contribution is also 
linear in T and has to be kept to be consistent. We em- 
phasize that the use of a simple rate-equation picture, 
which neglects this exchange field and takes into account 
spin-dependent tunneling only, would be fundamentally 
insufficient. 

The type of exchange interaction that emerges in our 
theory has been discussed in the literature in the context 
of Kondo physics for magnetic impurities in (normal) 
metalsiSl With the help of a Schrieffer- Wolff transfor- 
mation, the Anderson Hamiltonian describing the mag- 
netic impurity can be transformed to an s-d model, in 
which the spin of the magnetic impurity is coupled to 
the conduction-band-electron spins of the metal. Since 
the metal is spin symmetric, there is no net exchange 
field, though. This contrasts with our model, which in- 
cludes a finite spin polarization of the leads. By integrat- 
ing out the lead electrons of the transformed Hamiltonian 
in the subspace of single dot occupancy, we recover the 
precise mathematical form of the exchange field as given 
in Eq. (j^ . 

The generalization of our theory to the case of real fer- 
romagnets with nonuniform density of states is straight- 
forward. The only modification is to keep the energy de- 
pendence 0fpr{uj) = [pr+{(jj) — pr-i(^)]/[Pr+{i^)+Pr-{l-^)] 

and Tr{Lij) in all formulas. In particular, these quantities 
have to be included in the energy integral for the ex- 
change field in Eq. (ESJ.^^'^^ 

We remark that the effective exchange field is not only 
responsible for the torque on an accumulated spin as dis- 
cussed above, but it also generates a spin splitting of the 
dot leveh'^^ Since this Zeeman-like splitting is propor- 
tional to r, it cannot be resolved in any transport signal 
in the weak-coupling regime considered in the present 
paper. This splitting gives rise to a correction of higher 
order in the coupling that has to be dropped in a con- 
sistent first-order-transport calculation. This contrasts 
with the limit of strong lead-dot coupling, in which any 
spin splitting modifies the emerging Kondo physics in a 
crucial manner 



reservoir r. 



du; tr{r,/+MG>H + r,/-HG<(c^)} 



(4.1) 



Here, the bold face indicates a 2 x 2 matrix struc- 
ture for the spin degree of freedom. The lesser and 



greater Green's functions G^(aj) 



G<(a;) G< (c.) 



r 



and the same for Fr but 



and G^(w) are defined as the Fourier transforms of 
G<,,it) = licliOMt)) and G>,,(i) = -j(c.(i)ct,(0)), 
respectively. The coupling matrices are defined as Fl = 
I PL e-'-^/^ 
PL 1 
with the substitutions L ^ R and ((> ^cf). The diago- 
nal matrix elements describe tunneling of spins along the 
z direction and are, therefore, independent of the leads' 
spin polarization since the magnetizations of the leads 
lie in the x — y plane. This contrasts with the off-diagonal 
matrix elements taking into account tunneling of spins in 
the x — y plane. 

In the weak-coupling limit, the current is dominated by 
the contributions of first order in the coupling V. Because 
r already explicitly appears in Eq. 1)4. I|l . the Green's 
functions have to be calculated to zeroth order in Fi^i 
The final result for the current can be written as a sum 
over two contributions, 

I, = /P -t- /S , (4.2) 

which explicitly contain either the occupation probabili- 
ties Poj^ii-Pd or the accumulated spin S, respectively. 



/+(£ + !/) -/-(e) 



Pi 



-.f-{e + U)Pd 



(4.3) 



r 



-Pr Er 



2(-e) 



[f-ie) + f+{e + U)]S-hr. (4.4) 



Eor nonmagnetic leads, only the term contributes. At 
finite spin polarization in the leads, the additional term 
if appears, but also the first contribution is affected 
via the modification of the occupation probabilities as 
given in Eq. 1)3. 4|l . 

In steady state the two currents /l and /r are con- 
nected by the conservation of charge on the dot, so the 
homogeneous current flowing through the system can be 
defined as / = /r = — /l. 



V. RESULTS 



IV. ELECTRIC CURRENT 

The current through the device can be expressed 
in terms of Keldysh Green's functions of electrons, as 
shown, e.g., in Refs. 24"3^ The starting point of our 
discussion is the relation for the current Ir coming from 



In the following we discuss results for the spin accu- 
mulation and its impact on transport for both the linear 
and nonlinear regime. To be specific, we choose sym- 
metric coupling Fl — Fr — F/2, equal spin polariza- 
tions Pl = Pr = a-iid a symmetrically applied bias 
Vn = -Vi^ = V/2. 
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A. Linear response 

In equilibrium, i.e., without any applied bias volt- 
age, y = 0, no current flows and the stationary solu- 
tion of the rate equations for the occupation probabili- 
ties (EH) is given by the Boltzmann distribution, ~ 
exp(— S^/fceT), and since in this case ((i5'/di)acc = 0, no 
spin accumulation occurs in the quantum dot, = Oi^ 

At nonzero bias voltage, a current flows through the 
system and spin is accumulated on the dot. To describe 
the dot state in the linear-response regime, eV <^ fceT, 
we expand the steady-state solution of Eqs. (|3.4(l and 
(|3.5|) up to linear order in V. We find that (due to 
symmetric coupling) the dot occupation probabilities 
(Pqj -Pi) Pd) have no corrections linear in V and are, thus, 
unchanged, independent of the size and direction of the 
leads magnetization. In contrast to the occupation prob- 
abilities, the source term for the spin accumulation has 
a nonzero contribution linear in V, 

vT —eV , ^ , , 

generating a finite spin polarization along hl — riR, i.e., 
along the y axis. This is consistent with the notion of an- 
gular momentum conservation when magnetic moments 
polarized along hl enter from the left than those along 
fiR leave the dot to the right. 

The damping term, {dS/dt)i-ci — 
-(r/fi) [/-(e) +f+{e + U)] 5, limits the magnitude of 
spin accumulation, but it does not affect its direction. 
The third term, however, {dS/dt)rot = S x (Sl + Br) 
with Br - nrPr/{2TTh) j' duj (^j^^ - ^) f {oj) and 

Bl + Bii = Bo cofi{(p/2) Bx, yields a precession of the 
spin with an angle 

a = — arctan ^BqTc cos (5-2) 

about hl +nR, i.e., the a;-axis. As a result, the accumu- 
lated spin acquires both y and z components as seen in 
Fig. 01 This precession acts back on the magnitude of the 
accumulated spin. The stationary solution of Eq. (|3.10(l 
yields \S\ = Tc\{dS/dt)^^J cos a, from which we infer 
that the precession is accompanied with a decrease of 
the spin accumulation. 

The rotation angle a is plotted in Fig. IHb) as a func- 
tion of the level position e. We see that a changes 
sign at e = ~U/2, which reflects a sign change of the 
exchange field at this point. The point e = —U/2 is 
special, as in this case the system bears particle-hole 
symmetry. The fast variation of a in this region gives 
rise to a rather sharp feature in the magnitude of ac- 
cumulated spin as a function of e, see Fig. EJc). The 
width of the emerging peak scales (for < U /2) with 
i:U[f~{e) + f^{e + U)]/[p cos((/)/2)], which can be even 



X 




FIG. 4: Spin dynamics in the linear-response regime. Spin 
accumulates along the y direction. The spin precesses due to 
the exchange field that is along the x direction. Therefore, 
the stationary solution of the average spin on the dot is tilted 
away from the y axis by an angle a, plotted in Fig. |3b). 



smaller than temperature since the factor with the Fermi 
functions involved can become numerically small. 




FIG. 5: (a) Linear conductance normalized by r/fceT as a 
function of the level position e for different angles <j). (b) Angle 
a enclosed by the accumulated spin and the y axis as defined 
in Fig. 21 (c) Derivation of the magnitude of the accumulated 
spin on the dot with respect to the source-drain voltage V . 
Further parameters are p = 0.9 and U = IQk^T. 



As pointed out above, in the regime under considera- 
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tion the occupation probabilities do not depend on the 
spin polarization of the leads. In particular, they are 
independent of the relative angle (j> of the leads' magne- 
tization. This means that the (jj dependence of the con- 
ductance will be completely determined by the product 
S hi, = — S'-nR, as can be seen from Eqs. (14.3(1 and 1(4. 4|l . 
It is the relative orientation of the accumulated spin and 
the drain (or source) that produces the dependence of 
the current, rather than the product nL-nR, as in the 
case of a single magnetic tunnel junction. In this way, 
the dependent hnear conductance G^™ = {dI/dV)\v=o 
directly reflects the accumulated spin. The effect of the 
exchange field B is seen from the analytic expression 



Gi'"(0) 



= 1 



P 



sin^((/)/2) 



l + (Bore)2cos2(0/2) ' 



(5.3) 



which is plotted in Fig.|H|for different values of the level 
position e. 
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FIG. 6: Normalized conductance as a function of the angle 
enclosed by the lead magnetization for different level positions 
and the parameters U = IQk^T and p = 0.9. 



B. Nonlinear response 

We now turn to the discussion of the non-linear re- 
sponse regime, eV > k-^T . In Fig. Cfa) we show the 
current / as a function of the bias voltage V for an an- 
tiparallel configuration of the leads' magnetizations and 
different values of the leads' spin polarization p. As it is 




eV/kgT 



For e > 0, the quantum dot is predominantly empty, 
and for £ + [/ < doubly occupied with a spin singlet. 
This is reflected by a small lifetime Tc of a dot spin. Due 
to the short dwell time of a single spin-polarized electron 
in the dot, the rotation angle is small and the normalized 
conductance as a function of the relative angle of the 
lead magnetizations shows a harmonic behavior, see, e.g., 
the curve for e = SfceT in Fig. 

For —U < e < 0, however, the dwell time is in- 
creased and the exchange field becomes important. It 
causes the above described spin precession, which de- 
creases the product S'-hl since the relative angle <(nL, S) 
is increased and the magnitude of S is reduced. Accord- 
ing to Eq. 1(4.4(1 . the spin precession, thus, makes the 
spin-valve effect less pronounced, leading to a value of 
the conductance that exceeds the expectations made by 
Slonczewsk:^ for a single magnetic tunnel junction. 

For parallel and antiparallel aligned lead magnetiza- 
tions, (j) — and = TT, the accumulated spin and the 
exchange field also get aligned. As a result, there is no 
spin precession, even though the exchange field is still 
present. For these cases, the (^-dependent conductance is 
not affected by the exchange field, see Fig. |S1 



FIG. 7: Current-voltage characteristics for antiparallel (a) 
and perpendicular aligned (b) lead magnetizations. Further 
parameters are Fl — Fr = r/2, Pl = Pr = p, e = lOksT, 
and U = aOfeT. 

well known from transport through quantum dots with 
nonmagnetic leads, the current increases in a steplike 
manner. At low bias voltage, the dot is always empty and 
transport is blocked. With increasing bias voltage, first 
single and then double occupancy of the dot is possible, 
which opens first one and then two transport channels. 
A finite spin polarization p leads to spin accumulation 
and, thus, to a reduction of transport. A reduction of 
transport with increasing p is also seen for noncoUinear 
magnetization. There is a qualitative difference, though. 
As can be seen in Fig.|3b), a very pronounced negative 
differential conductance evolves out of the middle plateau 
as the spin polarization p is increased. 

Before explaining the origin of this negative differential 
conductance, we address the height of the plateau at large 
bias voltages first. Away from the steps, all appearing 
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Fermi functions are either or 1. From Eq. H4.2|l we get 
/ = {eT/2h) [1 — pS ■ (hl — H-r)]- When the exchange 
field is ignored, which can be done here since the spin 
dwcU time is shortened by the possibihty of forming spin 
singlets, then the accumulated spin is S — p{nL — nR)/4 
from which we get 



2n 



2 ■ 2 V 

p sm — 
^ 2 



(5.4) 



The suppression of transport due to the spin polarization 
p of the leads is comparable with the case of a single- 
tunnel junctions, when charging effects are of no impor- 
tance. 

At intermediate bias voltages, the dot can be empty or 
singly occupied, but double occupation is forbidden. To 
understand the negative differential conductance we first 
neglect the exchange field and then, in a second step, an- 
alyze how the exchange field modifies the picture. Since 
double occupation of the dot is prohibited, all electrons 
entering the dot through the left barrier find an empty 
dot. This is consistent with the fact that the current 
I = (eT/h)PQ explicitly depends only on the probability 
to find the dot empty. The spin accumulation affects the 
current only indirectly via modifying Pq. In the absence 
of the exchange field the relation between them is 



S=p 



— PoHL Z Hr 

i R ^ 



1 



1 + 2rL/rr 



with 



rR/rL + 2 i-p2 



(5.5) 
(5.6) 



(At this point we explicitly allow for different tunnel- 
ing strengths.) As a consequence, with increasing (f) and 
p the probability Po decreases, i.e., the dot is most of 
the time occupied with one electron, and the accumu- 
lated spin tends to align along — Hr, antiparallel to the 
drain electrode. The electron is, thus, trapped in the 
dot since tunneling out to the drain electrode is sup- 
pressed by a low density of states for the given spin direc- 
tion. No second electron can enter the dot because of the 
Coulomb interaction, and transport is, thus, blocked. So 
this mechanism is a type of spin blockade but with a dif- 
ferent physical origin compared to the systems described 
in literaturei^^iS^ 

The suppression factor coming from Eq. H5.6|l defines 
the local minimum of the current in Fig. El^b)- At this 
point, the relevant exchange field component generated 
by the coupling to the left lead vanishes, so that spin pre- 
cession becomes insignificant, see Fig. |^a). Away from 
this point spin precession sets in as illustrated in Fig. |S1 
The spin rotates about iiL and the electron can now more 
easily leave the dot via the drain electrode. 

Similar to the linear-response regime, the rotation 
caused by the exchange field weakens the spin-valve ef- 
fect. The two mechanisms, the spin blockade, and the 
spin precession have opposite effect on the spin-valve be- 
havior. Together with the fact that the strength of the 




FIG. 8: For electrons polarized antiparallel to the drain lead, 
the influence of the effective field generated by the source 
lead is dominating. By rotating the spins, the spin blockade 
is lifted and therefore the conductance recovers. 



exchange field varies as a function of the level position 
relative to the Fermi level, we now fully understand the 
origin of the nonlinear conductance. To illustrate this 
further, we plot in Fig.|5fb) the current which we obtain 
when we drop the contribution, Eq. (|3.8(l , of the exchange 
field to the spin dynamics by hand and compare it with 
the full result. In the absence of the exchange field, a 
wide plateau is recovered, whose height is governed by 
Eq. (|5.t)|l . The peak at the left end of the plateau indi- 
cates that, once the dot level is close to the Fermi level 
of the source electrode, the spin blockade is relaxed since 
the dot electrons have the possibility to leave to the left 
side. 

We now estimate under which circumstances the neg- 
ative differential conductance will be visible. Since the 
spin-blockade mechanism is crucial for the negative dif- 
ferential conductance to form, the relevant criterion can 
be obtained from Eq. H5.6I) as 



rR/rL + 2 1-p 



P ■ 2'P 

— 2 



(5.7) 



i.e., for a typical value (/) « 7r/2 and symmetric coupling 
Fl = Fr we need at least a spin polarization p 0.77, 
in agreement with Fig. [71 An asymmetry in the cou- 
pling strength Fl > Fr helps us to somewhat reduce the 
number of the required spin polarization. 

The effect of the spin blockade on the (/)-dependent 
current is depicted in Fig.^| We choose the bias voltage 
according to eV/2 — e + U/2, such that the influence of 
the exchange field is absent. For p — 0.5 still a sin^ ^ 
dependence can be recognized. For higher polarizations 
the conductivity drops faster and stays nearly constant 
at its minimal value due to spin blockade. This is just the 
opposite behavior than predicted for the linear-response 
regime as seen in Fig. 

We close this section with the remark that while we 
plotted only results for the case e > 0, in the opposite 
case e < the current-voltage characteristics is qualita- 
tively the same. 
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eV / k„T 



FIG. 9: Panel (a) The absolute value of the effective exchange 
field contributions from the left and right leads. Panel (b) the 
current voltage dependence, with and without the influence 
of the exchange field. For both plots the parameters (j) = n/2, 
Vl = Vr = r/2, e = fOfesT, U = SOfesT, and p = 0.95 were 
chosen. 



VI. EXTERNAL MAGNETIC FIELD 

As we discussed in the previous section, the trans- 
port properties through the quantum dot crucially de- 
pend on the magnitude and direction of the average 
spin in the dot, which is determined by the interplay 
of spin-dependent tunneling and the effective exchange 
field. The application of an additional, external mag- 
netic field Bcxt opens the possibility to manipulate the 
dot spin and, thus, the transport behavior in a direct 
way. Since the external magnetic field is independent of 
all the system parameters considered so far, the parame- 
ter space for experiments is enlarged by another dimen- 
sion. Magnetic fields have been used either to split the 
energy levels for the different spin components'^ or to 
rotate an accumulated spinii^ A combination of both, a 
static field splitting the energy levels plus an oscillating 
field to rotate the spin, defines an ESR scheme that al- 
lows for the study of single-spin dynamics in quantum 
dots in the absence of ferromagnetic leads, as proposed 
in Ref.Ei 




" ^..^ p=0.95 , - " 

aI 1 \ 1 A_ ^ 1 ^ 

TZil ^ 371/2 271 

FIG. fO: Angular dependence of the conductance with an 
applied voltage oi V = e -\- U/2, i.e., the voltage generating 
the smallest influence of the exchange field. Further plot pa- 
rameters are Fl = Fr = F/2, Pl = Pr = p, e = IQkBT, and 



The way a static external magnetic field affects the dot 
qualitatively depends on the ratio of the induced Zeeman 
energy and the intrinsic line width F of the dot energy 
levels. 



A. Large Fields 

At large magnetic fields, Bcxt S> F (to keep the for- 
mulas transparent, we include the gyromagnetic factor 
in the definition of i^cxt), all the complex spin dynam- 
ics discussed so far is absent. It is natural in this case 
to quantize the dot spin along the direction of Bcxt- 
The Zeeman splitting term i?oxt = ~ i'^ Li- 
ouville equation (|3.2|l ensures that for weak coupling, 
the density matrix always remains diagonal, as can be 
seen from expanding Eq. (|3.2I) to lowest order in F to 
obtain = (e| — £i)P\ ■ As a consequence, the aver- 
age spin on the dot is always polarized along the di- 
rection of the external magnetic field, and the exchange 
field does not play any role. In a physical picture, the 
spin components perpendicular to Boxt are averaged out 
due to a fast precession movement. Transport effectively 
takes place through two energetically separated levels, 
where the coupling strength to the left and right leads 
depends on the angles 6^ — <(Boxt,nr) between the 
lead magnetizations and the external field via the rela- 
tions Fr| ~ Vr{l+Pr COS Or) and Trl = T r(\ — Pr COS Or)- 

Thereby we assumed that the applied magnetic field is 
not strong enough to tilt the lead magnetization direc- 
tion. It is worth to point out that the spin polarization 
and the trigonometric function of the relative angle Or ap- 
pear always as a product, so that this system is equivalent 
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to the colhnear casesiSiiS with decreased magnetization 
of the leads. 



B. Small Fields 

The situation becomes more interesting in the limit of 
small fields, i?cxt S T. To describe this case, we per- 
form the same perturbation expansion in F as before, 
but count -Bcxt as being first order in T. This introduces 
in Eq. H3.2|l the (first-order) term (£| — ei)Pj, while the 
kernels S are then not affected by the external magnetic 
field, because they are already of first order in the cou- 
pling strength T. 

It turns out that, in the considered limit, the only 
change introduced by the external magnetic field is an 
additional term S x Bext in the Bloch-like equation 



dS_ 



dS 

dt 



+^ +SxB 



dS 

dt 



(6.1) 



with B = Sl -|- Br -I- Bcxt- With the external magnetic 
field as a new parameter all the effects of spin blockade 
and the precession discussed so far can be enhanced or 
suppressed and made independent of the lead alignment. 

An external field is applied along ex = (nL + nR)/|nL-l- 
ArI leads to a suppression of the spin- valve effect in the 
linear-response regime similarly as due to exchange field. 
In fact, according to Eq. H6.1|l the magnitude of the rel- 
evant field is just the sum of the exchange and the ex- 
ternal fields. The linear conductance as a function of 
Bext is displayed in Fig. In the absence of the effec- 
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FIG. 11: Linear conductance as a function of an applied 
external magnetic field in x direction. The parameters are 
Fl = Fr = F/2, PL = PR = p, e = 0, [/ = WksT, and 

(j> = Tv/2. 



five exchange field the conductance would be symmetric 
around Bcxt = 0, and any value of the field would increase 
the conductance as discussed in Sec. IV Al The exchange 
field yields a shift of the position of the minimum. This 
opens the possibility to determine the exchange field ex- 
perimentally by changing the external magnetic field and 
keeping all other parameters fixed. From the generaliza- 
tion of Eq. (|5.3|) to the present case it becomes apparent 



that the width of the minimum is given by the inverse of 
the dwell time Tc in the dot, since during this time the 
dot spin is influenced by the fields. 

We remark that the effect of external-field components 
perpendicular to the x direction is somewhat more com- 
plicated, even in the linear-response regime. The reason 
is that the induced precession will generate some finite 
spin accumulation along the x direction, which intro- 
duces a linear correction in V of the charge occupation 
probabilities via Eq. H3.4|) . This, in turn, acts back on 
the accumulated spin due to Eq. ()3.6|) . 



VII. INCLUSION OF SPIN RELAXATION 

We emphasize that the only source of spin relaxation 
considered in our analysis is the tunnel coupling of the 
dot states to the left and right leads, i.e., we require 
that the spin-flip relaxation time due to other mecha- 
nisms such as coupling to nuclear spinS)^ or spin-orbit 
coupling2& is large enough not to destroy spin accumu- 
lation. If we assume a spin relaxation rate T2 equal to 
the ensemble spin dephasing time T2 ~ 100 ns for n- 
doped GaAs bulk material-- as lower boundary, it is still 
justifiable to neglect spin relaxation effects for interface 
resistances of the order of IM^l.'^- But recent measure- 
ments with semiconducting quantum dots indicate even 
much longer spin relaxation times. 

Nevertheless, (extrinsic) spin relaxation can be the 
most limiting factor for various spin-accumulation exper- 
iments in systems of larger size.i'^''*^' To include extrinsic 
spin relaxation in a simplified phenomenological way, we 
add an isotropic relaxation term — S/roxt in Eq- (|3.7|l . 
which leads to 



dS_ 

Itt 



rol 



1 

''ext 



(7.1) 



Now, two characteristic time scales have to be distin- 
guished: the average time an electron stays in the dot, 
Tc, and the meanlife time of a dot spin, Tg = (l/'''c + 

1/Tcxt)^^ 

By solving the system of master equations for the 
linear-response regime, we see that we have to replace 
the factor with Tg everywhere in the previous discus- 
sions, which matches the intuition that Tg is the rele- 
vant time scale for spin related effects. From Eq. (|3.10() 
we see that the extrinsic relaxation processes reduce the 
magnitude of the spin accumulation by a factor Ts/tc. 
The observation that the accumulated spin on the dot 
and, therefore, all spin related effects in transport vanish 
linearly with Ts/tc (and not exponentially) implies that 
magnetoresistance effects can still be measured, even if 
the mean dwell time of the electrons in the dot exceeds 
the spinlife time significantly, which was experimentally 
confirmed by Ref. l42l The reason for the linear depen- 
dence is the following. Even when the spins relax much 
faster than the mean dwell time, a fraction of the spins 
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leaves the dot already before being relaxed. The proba- 
bility of detecting a spin with the original alignment in 
a transport measurement is given by the time integral 
^spin = /q dtexp(—t/Ts)/Tc — Ts/tc, which is linear de- 
spite the exponential decay of the accumulated spin on 
the dot. 



VIII. RELATION TO SPIN-MIXING 
CONDUCTANCE 

Transport through multiterminal devices with non- 
collinear lead magnetizations has been discussed recently 
by using the language of a complex spin-mixing conduc- 
tance G\i within a circuit theory.-^ The mixing conduc- 
tance G-\i contains the information about the transport 
of spins oriented perpendicular to the magnetization of 
the ferromagnetic leads. Using this approach, the con- 
ductance of a two-terminal device with noncoUinear mag- 
netization of the leads can be described in linear response 
by the general relation 



G{9) = 



G 



tan2(6'/2) 



tan2(6'/2) + \Tj\^/Re7j 



(8.1) 



where r] = 2G-\i/G is the complex spin-mixing conduc- 
tance normalized to the sum G of the spin-up and spin- 
down conductance of a single contact. This formula, 
Eq. (|8.1|l . is quite general but requires the knowledge 
of 77. Hybrid systems of ferromagnets with normal met- 
als in diffusive and ballistic regimeS^ or even Luttinger 
liquids^^*^ have been described. For system with tunnel 
contacts Re 77 = 1 and Eq. (|5.3(l can be transformed into 
the form Eq. (|8.1|l by identifying 



2ImG 



Ti 



G 



BnT. 



(8.2) 



Our approach, thus, provides a theory for the interac- 
tion effects on the spin-mixing conductance in the linear- 
response regime. The imaginary part Im 77 is proportional 
to the local effective field experienced by the dot spins, 
because this field is responsible for the spin precession 
which yields a coherent mixing of the two spin channels. 



IX. CONCLUSIONS 

In the presented work the influence of charging interac- 
tion effects on the tunnel magnetoresistance for systems 
with noncoUinearly magnetized leads has been discussed. 
We examined a single-level quantum connected to ferro- 
magnetic leads as a model system and found that charge- 
and spin-related properties do not just coexist but their 
interplay gives rise to new effects detectable in transport 
measurements. 

The most important effects discussed in this paper are 
related to the existence of an effective exchange magnetic 
field Bl + -Br, which leads to a precession movement of 



the accumulated spin on the dot. This precession weak- 
ens the spin- valve effect by reducing the amount of accu- 
mulated spin and changing its orientation relative to the 
leads. As a result, the current flowing through the system 
is increased due to the exchange field. Also the functional 
form of the angular-dependent linear conductance is af- 
fected, and the cos-dependence valid for a single magnetic 
tunnel junction is modified due to the precession of the 
accumulated spin. 

In nonlinear response, the accumulated spin on the dot 
tends to align antiparallel to the drain lead magnetiza- 
tion. Such a spin blockade causes a suppression of trans- 
port for the bias-voltage regime in which the dot can 
be empty or singly occupied. The spin precession due 
to the coupling to ferromagnetic leads countersteers the 
spin blockade by tilting the spin accumulation direction. 
Since the strength of the exchange field is a nonmonotonic 
function of the applied voltage, a regime with negative 
differential conductance exists. 

Attaching ferromagnetic leads to quantum dots is quite 
challenging. One promising approache is to contact an ul- 
trasmall aluminum nanoparticle, which serves as a quan- 
tum dot, to ferromagnetic metallic electrodes. In this 
way, quantum dots with one magnetic and one nonmag- 
netic electrode have already been fabricated. An al- 
ternative idea is to use carbon nanotubes as quantum 
dots and to place them on ferromagnetic contactsji^* But 
also other systems are conceivable, such as self-assembled 
InAs quantum dots in GaAs LEDs with (Ga,Mn)As as 
ferromagnetic electrode, "'^^ or STM setups with a natural 
or artificial impurity on a ferromagnetic substrate con- 
tacted by a ferromagnetic STM tip4& 
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APPENDIX A: DIAGRAMMATIC TECHNIQUE 

In this appendix we present some technical details of 
the derivation of the master equations and the transport 
current. The approach is based on the diagrammatic 
technique developed in Ref. |2^ designed for the applica- 
tion in quantum dot systems. The actual calculation is 
to some degree similar to transport through Aharonov- 
Bohm interferometers with quantum dotsi^ 
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Transition rates E 



In the Liouvillc equation (|3.2I) the transition rates ap- 
pear as irreducible self-energies S^J . These self-energy 
parts are represented as block diagrams enclosed by a 
Keldysh contour. Examples of first-order diagrams are 
shown in Fig. 1121 The real-time axis runs horizontally 
from the left to the right while the upper (lower) line 
represents the forward (backward) propagator of the dot 
on a Keldysh time contour. The rules to calculate the E 



0- 



^dd 



T 



(i)..^^' + (T) -^^y.^ ^) 

Q-^ i ^ i 




d d 

d d'^ — 



d 



(T)"~l>'' (i) + (i)- (T) 
T I I ] 



over all intermediate lead states (r = R/L, cr = +/— ). 
Dependent on the time ordering, the transition rate is 
proportional to the electron /+ or the hole distribution 
function /~ — 1 — /+, while the relative phases can be 
extracted from Fig. |3| 

5. ) The diagram gets a prefactor of (— 1)^(— 1)"^, where 
h is number of internal vertices on the backward propa- 
gator, and c the number of crossings of tunnel lines. 

6. ) Integrate over all energies of the tunneling lines. 
In the sequential-tunneling regime, only one tunnel line 

per diagram is involved. Therefore, only one frequency 
integral appears which can be calculated trivially by us- 
ing Cauchy's formula. 



2. Green's functions 

For the calculation of the transport current, we need 
the Fourier transform of the Keldysh Green's functions, 

G>,,{t) = -*(c,(04,(0)), (A3) 
G<At) = +^(4'(0)c,(^)). (A4) 
The Green's function is an average of one creation and 



FIG. 12: Examples for the diagrams representing the gener- 
alized transition rates E. 

areii^ 

1. ) Draw all topologically different diagrams with tun- 
neling lines connecting vertices on either the same or op- 
posite propagators. Assign to the four corners and all 
propagators the states x = 0, t, i, d and the correspond- 
ing energies e^, as well as an energy oj to each tunneling 
line. 

2. ) For each time interval on the real axis confined by 
two neighboring vertices, assign the resolvent l/(Ai? -|- 
10+), were is the energy difference between left and 
right going tunnel lines and propagators. 

3. ) For each vertex connecting a double-occupied state 
d to the up state the diagram acquires a factor (—1). 

4. ) Assign to each tunneling line the factor 7^^, (cj) 
[7^7^/ {lo)] when the tunneling line is going backward (for- 
ward) with respect to the Keldysh contour. Its form is 

7±(c.) = ^/.^M- (Al) 

r,x=± 

-rL-/±e'*/2 - rR_/±e-**/2 ) (A2) 

Here, a and a' are the spins of the electron that leaves 
and enters the vertices connected by the line. The factors 
7^(0;) come from the contraction of two lead operators 
in the tunnel Hamiltonian and resemble the transition 
rates predicted by Fermis golden rule, including a sum 




a 



FIG. 13: Examples for the diagrams representing Green's 
functions. 

one annihilation operator. These operators appear as 
external vertices (open circles) in the diagrams. They 
can be constructed following the rules: 

1.) Draw the bare contour, the external vertices, and 
the virtual tunnel line in the proper orientation as seen 
in Fig. El The indices of the Green's function indicate 
the spin of the created and annihilated electrons at the 
external vertices. If higher-order contributions in F are 
to be calculated, add internal vertices according to the 
desired order and connect them by tunnel lines. Multi- 
ply the diagram with the matrix element of the density 
matrix corresponding to the inital state. Assign to each 
propagator with state x the corresponding energy e-^ and 
to each tunnel line and the virtual line an energy uj. 
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2. ) Apply the same rule as no. (2) for calculating S. 

3. ) Apply the same rule as no. (3) for calculating E. 

4. ) Assign each tunneling line the factor {—ly^^^, (uj), 
where v is the number of external vertices on the Kcldysh 
contour enclosed by the end points of the tunnel line. The 
line connecting the external vertices does not contribute. 

5. ) Apply the same rule as no. (5) for calculating S. 

6. ) Integrate over all energies except the energy as- 
signed to the virtual tunnel line. 



For the lesser function G^^, , the direction of the ex- 
ternal tunnel lines is reversed, and an additional global 
minus sign is involved in the definition. 

In this paper we assume weak dot-lead coupling, i.e., 
we need only the zeroth-order Green's functions. In this 
limit, we get G^^{llj) — —2TTiPg6{u! — £—U) — 2TriPod{LLj — 
e), G<„{cj) = +2TTiP„Siuj-e) + 2TTiPdS{uj-e-U) and for 
CT 7^ a we find G>^{uj) = 2mP§5{uj -e-U), G<s{lu) = 
2mP§5{uj-e). 
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